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Abstract 

The two-point autocorrelation function of ultra-high energy cosmic ray (UHECR) 
arrival directions has a broad maximum around 25 degrees, combining the data 
with energies above 4 x 10^® eV (in the HiRes energy scale) of the HiRes stereo, 
AGASA, Yakutsk and SUGAR experiments. This signal is not or only marginally 
present analyzing events of a single experiment, but becomes significant when data 
from several experiments are added. Both the energy dependence of the signal and 
its angular scale might be interpreted as first signatures of the large-scale structure 
of UHECR sources and of intervening magnetic fields. 

RAGS: 98.70.Sa 


1 Introduction 


The sources of ultra-high energy cosmic rays (UHECR) are despite of more 
than 40 years of research still unknown. Main obstacle for doing charged parti¬ 
cle astronomy are deflections of the primaries in the Galactic and extragalactic 
magnetic helds. While the magnitude and the structure of extragalactic mag¬ 
netic fields are to a large extent unknown, already deflections in the Galactic 
magnetic field alone are large enough to prevent UHECR astronomy if the 
primaries are heavy nuclei [1,2]. Assuming optimistically that the primaries 
are protons, typical deflections in the Galactic magnetic field are around five 
degrees in most part of the sky at E = 4 x 10^® eV [2]. Therefore, it might 
be possible to perform charged particle astronomy, if moreover deflections in 
extragalactic magnetic fields are sufficiently small. 

This scenario can be divided in two quite different sub-cases: In the hrst one, 
a small number of point sources results in small-scale clusters of arrival direc¬ 
tions around or near the true source positions. Accumulating enough events. 


Preprint submitted to Elsevier Science 


20 December 2005 




the identification of sources will become possible using e.g. correlation studies. 
Various studies have been pursued in this direction [3,4], In the second sub¬ 
case, the number of point sources is too large to receive two or more UHECRs 
from the same source with the present statistics. However, the sky density of 
UHECRs reflects the large-scale structure of the sources and, possibly, of the 
intervening magnetic helds [5,6]. Therefore, the measured UHECR distribu¬ 
tion is anisotropic and over-/underdense regions exist that reflect the angular 
size of up-to 15-20 degrees of typical structures in the galaxy distribution. 
Obviously, Nature might have chosen a mixture of these two extreme possibil¬ 
ities: The vast majority of UHECR sources might produce only singlet events, 
while a subclass of sources with extreme luminosity might be detectable as 
point sources via small-scale clustering studies. Furthermore, point sources 
might be easier to identify at the highest energies, if the number density of 
sources decreases with the maximal energy Umax to which they can accelerate 
as argued in Ref. [7]. 

In this work, we study the arrival direction distributions of the UHECRs, 
putting emphasis in contrast to most earlier studies on intermediate angular 
scales. Since these two-dimensional distributions average three-dimensional 
structures (with typical scale L) over the mean free path I of UHECRs, no 
anisotropies reflecting the large-scale structure of sources are expected for 
I ^ L. To obtain an optimal compromise between the number of events used, 
the mean free path I of UHECRs and deflections in magnetic helds, it is impor¬ 
tant to use a consistent energy scale when combining different experiments. 
Therefore, we discuss hrst in Sec. 2 the publicly available UHECR data and 
how we rescale the data of different experiments to a common energy scale. In 
Sec. 3, we analyze then the arrival direction distributions of the used UHECRs 
and discuss in Sec. 4 our results, before we summarize in Sec. 5. 


2 UHECR data sets and their energy scale 


The available data set of UHECR events consists of 

(1) The publicly available AG AS A data set until May 2000 from Ref. [8], 
consisting of 57 events with E > 4 x 10^® eV and zenith angle 'd < 45°. 
The exposure of the AGASA experiment is 4.0 x 10^® nV s sr [8]. 

(2) The Yakutsk data as presented at the ICRC 2005 [9]: 34 events with 
energy E > 4 x 10^® eV and zenith angle '& < 60°. The corresponding 
exposure is 1.8 x 10^® m^ s sr [9]. 

(3) The SUGAR data with energy above E > 1 x 10^® eV and zenith angle 
■d < 70° from Ref. [10]. We follow Ref. [11] and use 4.0 x 10^® m^ s 
sr as approximate exposure for i? = 55°. From this value, we rescale 
the exposure approximately to 3.0 x 10^® m^ s sr for -d < 45° and to 


2 


5.3 X 10^® s sr for -d < 70°, respectively. 

(4) The HiRes stereo data set: We deduce the arrival directions of the events 
with energy R' > 1 x 10^® eV from Ref. [12]. Then we divide the data set 
into events with energy in the range E = (1 — 2) x 10^®, E = (2 — 4) x 10^® 
and E > 4 X 10^® eV using Ref. [13]. Its total exposure is 4.6 x 10^® s 
sr [12]. 

(5) From the Volcano Ranch, Haverah Park, Flye’s Eye experiments no de¬ 
tailed informations are available about their events. Therefore, we can 
use only the events with E > 10^° eV for which the arrival directions 
are given in Ref. [14]: four events from Haverah Park, and one both from 
Volcano Ranch and Flye’s Eye. 

We use as angular acceptance ri{6) of a ground array experiment at geograph¬ 
ical latitude b that observes showers with maximal zenith angle 'dmax (cf., e.g.. 
Ref. [15]) 

Q^max 

ri{6) (X j da cos('d) oc [cos(5) cos(5) sin(amax) + ttmax sin(5) sin((5)] , (1) 

0 


where 


arccos(^) 

for 


TT 

for 

e<-i, 

0 

for 

e>i, 


( 2 ) 


and 


cos(^max) " sin(5) sin((5) 
cos(6) cos((5) 


( 3 ) 


In contrast to ground arrays, the exposure of Hires in right ascension R.A. 
is non-uniform. Therefore, we sample the exposure of HiRes from Fig. 2 of 
Ref. [12]. 

The absolute energy scale of each experiment has a rather large uncertainty. 
To reproduce correctly spectral features like the dip, the energies E given by 
the experiments have to be shifted to new energies E' [16]. First, we assume 
that the normalization of the HiRes stereo spectrum is consistent with the one 
of HiRes in monocular mode, following Ref. [17]. In Ref. [19], we compared 
the SUGAR data given in Ref. [10] to the energy spectrum measured by the 
AG AS A and HiRes experiments. We found that rescaling the SUGAR energies 
calculated with the Hillas prescription by 15% downwards, E' = UHi Has /Uld, 
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Fig. 1. Skymap of the UHECR arrival directions of events with rescaled energy 
E' > 4 X 10^® eV in equatorial coordinates; magenta crosses-30 Agasa (A) events 
with E > 5.2 X 10^® eV, red circles-27 HiRes (H) events with E > A x 10^® eV, 
black stars-13 Yakutsk (Y) events with E > 6 X 10^® eV, blue boxes-31 Sugar 
(S) events with E > 6 x 10^® eV, magenta crosses-4 Haverah Park (HP) events 
with E > 10^*^ eV, red triangle-one Flye’s Eye (FY) event with E > 10^*^ eV, blue 
triangle-Volcano Ranch (VR) event with E > 10^*^ eV. 

makes their data consistent with the ones from AGASA. The same is true for 
the HiRes spectrum if the energy is rescaled up-wards by 30%. In contrast to 
Ref. [19], we hx the energy scale by the HiRes mono data in the present work. 
Therefore, we have to shift the AGASA data by 30% downwards, and the 
SUGAR data by 50% downwards. According to Ref. [9], the Yakutsk energy 
scale is systematically 15-20% above the AGASA energy scale. Thus, in order 
to match the Yakutsk data to the HiRes energy scale we rescale all energies 
of UHEGR events of Ref. [9] by 50% downwards. 

In Fig. 1, we show a skymap in equatorial coordinates of the arrival directions 
of the UHEGR used in the analysis below. An inspection by eye indicates 
an overdense region around and south the AGASA triplet as well as several 
underdense regions or voids. In Fig. 2, we show for comparison a skymap 
with events from those two experiments (Hires and SUGAR) which published 
also arrival directions below U' = 4 x 10^® eV. Again an inspection by eye 
indicates that the addition of low-energy data appears to make the sky map 
more isotropic. In the next section, we perform a statistical analysis to deduce 
the typical angular scales of excess correlations visible for E' = 4 x 10^® eV. 
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Fig. 2. Skymap of the UHECR arrival directions in equatorial coordinates from 
Hires and SUGAR in three different energy bins, E = (1 — 2) x 10^® (magenta, 
small circles), E = (2 — 4) x 10^® (blue, medium box) and E > 4 x 10^® eV (red, 
large stars). 

3 Autocorrelation analysis 


We use as our statistical estimator for possible deviations from an isotropic 
distribution of arrival directions the angular two-point auto-correlation func¬ 
tion w. We define w as function of the angular scale 6 as 

N i-l 
i=l j=l 


where 0 is the step function, N the number of CRs considered and 6ij = 
acos (cos pi cos pj + sin pi sin pj cos(0j — 0j)) is the angular distance between 
the two cosmic rays i and j with coordinates (0, p) on the sphere. Having 
performed a large sample of Monte Carlo simulations, we call the (formal) 
chance probability P{5) to observe a larger value of the autocorrelation func¬ 
tion w{5) the fraction of simulations with w > w*, where w* is the observed 
value. We would like to warn the reader at this point that we have not fixed 
a priori our search and cut criteria. Thus the obtained probabilities are only 
indicative. But they can be used in particular to compare for different data 
sets the relative likelihood to observe the signal as chance fluctuation. 
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8 [degrees] 

Fig. 3. Chance probability P{6) to observe a larger value of the autocorrelation 
function as function of the angular scale 5 for different combinations of experimental 
data; label of experiments as in Fig. 1. 

In Fig. 3, we show the chance probability P{S) as function of the angular scale 
6 for different combinations of experimental data. The chance probability P{S) 
shows already a 2a minimum around 20-30 degrees using only the 27 events 
of the HiRes experiments with E' > 4 x 10^® eV. Adding more data, the 
signal around 6 = 25° becomes stronger, increasing from ~ 2a for 27 events 
to ~ 3.5cr for 107 events. It is comforting that the position of the minimum 
of P(S) is quite stable adding more data and every additional experimental 
dataset contributes to the signal. Moreover, autocorrelations at scales smaller 
than 25° become more signihcant increasing the dataset. However, we warn 
the reader again that we have not constrained ourselves a priori to search 
for autocorrelations at 5 = 20° and hence the numerical values of the chance 
probability found should not be taken at face value. 

To understand better how the search at arbitrary angular scales influences the 
signihcance of our signal we have calculated the penalty factor ^ for the scan 
of P(<5) over 6. The penalty factor increases for increasing resolution AS of the 
angular scale 6, but reaches an asymptotic value for AS —> 0. The numerical 
value of the penalty factor found by us in the limit Ad 0 varies between 6 
for the HiRes data set alone and 30 for the combination of all data. Since the 
energy cut we use is determined by the one chosen in Ref. [13], no additional 

^ For a discussion of the use of penalty factors see e.g Ref. [20]. 
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Fig. 4. Chance probability P{6) to observe a larger value of the autocorrelation 
function as function of the angular scale 6 for different cuts of the rescaled energy 
E': black E' > 1 x 10^® eV, blue E' > 2 x 10^® eV and red line > 4 x 10^® eV. 


penalty factor for the energy has to be included. We conclude therefore that 
the true probability to observe a larger autocorrelation signal by chance is 
P^3 X 10 ^ for the complete data set. 

In Fig. 4, we show again the chance probability P(^) as function of the angular 
scale 5, but now for different values of the minimal event energy included in 
the analysis. Again we can use only Hires and SUGAR data, since for the 
other experiments no arrival directions for events below E' = 4 x 10^® eV are 
published. Moreover, the energy bin size we use is dictated by the one chosen in 
Ref. [13]. The addition of data with energy below E' ^ 4x 10^® eV reduces the 
signihcance of the minimum of T’(5). This could have various reasons: First, 
the interaction length of the UHECR primaries can increase with decreasing 
energy, as in the case of protons or nuclei. Then, the projection on the two- 
dimensional skymap averages out more and more three-dimensional structures. 
Second, deflections in the Galactic and extragalactic magnetic fields destroy 
for lower energies more and more correlations. Note that the autocorrelation 
signal would not disappear lowering the energy threshold if it would be caused 
solely by an incorrect combination of the exposure of different experiments. 
In Fig. 4, there is also a prominent minimum of P(^) around 120° visible. 
This angular scale corresponds to the distance between overdense spots and 
is always for other datasets in Fig. 3 visible, but much less signihcant. 
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Fig. 5. Chance probability P{6) to observe a larger value of the autocorrelation 
function as function of the angular scale 5 combining the AGASA, HiRes, Yakutsk 
and SUGAR experiments according to the number of observed events (dashed lines 
labeled Nev) and according to the exposure of the individual experiments (solid 
lines labeled Exp). Data of SUGAR zenith angles up to 45 degrees and up to 70 
degrees are shown with thin and thick solid lines, respectively. 


In order to combine the various experiments in a single data set as shown in 
Figs. 3 and 4, we have estimated the relative exposure of each experiment by 
the number of observed events above a chosen energy threshold. We can check 
this procedure for those experiments for which we can cut the data set on the 
same energy, i.e. for AGASA, HiRes, Yakutsk and SUGAR. For this check, we 
hx the total number of events to the observed one, which is 92 events using 
zenith angle'd < 45° and 101 using 'd < 70° for SUGAR. Then we choose the 
event number of each experiment according to its exposure. For instance, for 
SUGAR with zenith angle "d < 45° we use 32 HiRes, 27 AGASA, 12 Yakutsk, 
and 21 SUGAR events. Note that the observed numbers of events are 27 
HiRes, 30 AGASA, 13 Yakutsk, and 22 SUGAR event, i.e. the differences 
are compatible with pure statistical (Poisson) fluctuations. The results of this 
comparison are presented in Fig. 5, combining the AGASA, HiRes, Yakutsk 
with SUGAR data with zenith angles up to 45 degrees and up to 70 degrees, 
respectively. The chance probability calculated with both recipes gives similar 
results. In one case the exposure method gives slightly stronger and in the 
other case slightly weaker results. 
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We have also checked how our result depends on the systematic energy dif¬ 
ference between the AGASA and HiRes stereo experiments. Up to now, we 
have assumed that the energies of the HiRes stereo data have no systematic 
shift relative to the HiRes mono data [17]. However, because the normaliza¬ 
tion of the HiRes stereo spectrum is still under discussion, we have checked 
as extreme possibility that the HiRes stereo energy scale is shifted upwards 
by 20%. In this case, the AGASA data should be shifted down only by 10%, 
while the Yakutsk and SUGAR data should be shifted by 26.5%. The total 
number of events in the sample will increase from 107 to 143. The resulting 
probability P{S) as function of the angle 6 is similar to one shown in Fig.3. 
Thus we conclude that our results are not strongly sensitive to the exact value 
of the systematic energy difference between the AGASA and the HiRes exper¬ 
iments. Note also that once the relative energy scales of the experiments are 
fixed, our results do not depends on value of the “true” energy scale, i.e. on 
the global shift of all experiments together up and down in energy. 


4 Discussion 


Our results, if conhrmed by future independent data sets, have several impor¬ 
tant consequences. 

Firstly, anisotropies on intermediate angular scales constrain the chemical 
composition of UHEGRs. Iron nuclei propagate in the Galactic magnetic field 
in a quasi-diffusive regime at U = 4 x 10^® eV and all correlations would be 
smeared out on scales as small as observed by us. Therefore, models with a 
dominating extragalactic iron component at the highest energies are disfavored 
by anisotropies on intermediate angular scales. 

Secondly, the probability that small-scale clusters are indeed from point 
sources will be reduced if the clusters are in regions with an higher UHEGR 
flux. For example, the AGASA triplet is located in an over-dense spot (cf. 
map in Fig. 1) and the probability to see a cluster in this region by chance 
is increased. In contrast, the observation of clusters in the ’’voids” of Fig. 1 
would be less likely by chance than in the case of an UHEGR flux without 
medium scale anisotropies. 

However, the most important consequence of our findings is the prediction 
that astronomy with UHEGRs is possible at the highest energies. The mini¬ 
mal energy required seems to be around E' = 4 x 10^® eV, because, as one 
can see from Figs. 2 and 4, at lower energies UHEGR arrive more and more 
isotropically. This trend is expected, because at lower energies both deflec¬ 
tions in magnetic helds and the average distance I from which UHEGRs can 
arrive increase. Since the two-dimensional skymap corresponds to averaging 
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all three-dimensional structures (with typical scale L) over the distance I, no 
anisotropies are expected for / ^ L. Thus, if the signal found in this analysis 
will be conhrmed it has to be related to the local large scale structure. 

Finally, we note that to check this signal an independent data set of order 
0(100) events with E' > Ax 10^® eV is required. This agrees roughly with the 
hnding of Ref. [6] that around 100 events are needed to reject the hypothesis 
that the UHECR sources trace the galaxy distribution. 


5 Summary 


We have found that the two-point autocorrelation function of UHECR arrival 
directions has a broad maximum around 25 degrees. Combining all publicly 
available data with energy E' > 4 x 10^® eV, the formal chance probability 
that a stronger autocorrelation is obtained from an isotropic distribution is 
around P ~ 1 x 10“^. Both the signal itself and the exact value of the chance 
probability have to be interpreted with care, since we have not hxed a priori 
our search and cut criteria. Taking into account a penalty factor of 30 appro¬ 
priate for our search at all angles 5 G [0 : 180°], we have estimated the true 
probability as P ~ 3 x 10“^. We have checked that the autocorrelation sig¬ 
nal disappears lowering the energy threshold, indicating that it is not caused 
solely by an incorrect combination of the exposure of different experiments. 
The autocorrelation signal found by us around 6 = 25° should be tested with 
future, independent data sets from HiRes, the Pierre Auger Observatory [21] 
and the Telescope Array [22]. If conhrmed, it constrains the UHECR primary 
type together with the magnitude of extragalactic magnetic helds and opens 
the door to astronomical studies with UHECRs. 
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